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We report 31 P-NMR and specific heat measurements on an iron (Fe)-based superconductor 
SrFe2(Aso.65Po.35)2 with T c =26 K, which have revealed the development of antiferromagnetic corre- 
lations in the normal state and the unconventional superconductivity(SC) with nodal gap dominated 
by the gapless low-lying quasiparticle excitations. The results are consistently argued with an un- 
conventional multiband SC state with the gap-size ratio of different bands being significantly large; 
the large full gaps in s±-wave state keep T c high, whereas a small gap with a nodal-structure causes 
gapless feature under magnetic field. The present results will develop an insight into the strong 
material dependence of SC-gap structure in Fe-based superconductors. 



Superconductivity (SC) in the 122 iron-pnictides 
MFe2As2 (M=Ba, Sr or Ca) can be induced by either 
electron- or hole-doping, an application of pressure or 
chemical substitution for M, Fe, and As. The substi- 
tution of K +1 for Ba +2 in (Bai_ 2 ,K a; )Fe2As2 (denoted 
as BaK122) suppresses antiferromagnetism (AFM) and 
structural transition through the hole-doping, and in- 
duces SC with the relatively high superconducting tran- 
sition T c =38 K at £=0.4 Extensive experimental 
works in this compound have revealed the fully-gapped 
s±-wave state [2-6]. The isovalent substitution of P for 
As in BaFe2(Asi_ x P x )2 (denoted as Bal22P) also ex- 
hibits the maximum of T c —30 K at £=0.33, providing 
a phase diagram similar to BaK122 without changing 
a carrier density [7[. Notably, a nodal-gap structure in 
the SC state has been reported by various experiments 
on Bal22P [7-12|. Relevant with this nodal-gap struc- 
ture, the 31 P-NMR study on Bal22P reported the 
gapless SC with a large fraction of residual quasiparticle 
density of states (RDOS) N ras at the Fermi level (Ep) 
with N TCS /N -0.34, where N is the DOS at E F in 
normal state. On the one hand, we note that a large 
fraction of N rcs /No —0.5 was reported even for the SC 
state of undoped SrFe2As2(Srl22) under highly hydro- 
static pressure [l^. It is therefore likely that the presence 
of RDOS is inherent in the 122 iron-pnictides in which 
a lattice contraction takes place by either the applica- 
tion of pressure or the chemical substitution of P for As. 
At the present, they are not clear which Fermi surfaces 
(FSs) are responsible for the existence of nodal-gap and 
why it exists there QO]. 

The replacement of either Sr or Ca for Ba site in 
A/122P compounds brings about further lattice contrac- 
tion. A phase diagram for SrFe2(Asi_ 2: P a ;)2 (Srl22P) 
resembles that for Bal22P[l8|, but a maximum T c de- 
creases to 26 K at £=0.35, and a maximum T c for 
CaFe 2 (Asi_ ;c P :E )2 (Cal22P) does to T c =15 K at x ~ 0.04 
in a narrow SC region[l9j. These results imply that the 



effect of lattice contraction may cause a reconstruction 
of FS topologies [lj}, causing the onset of SC to be un- 
favorable, but it provides us with a good opportunity to 
investigate an intimate change in electronic state and SC 
gap structure systematically in the M 122 iron-pnictides. 

In this paper, we present 31 P-NMR results on single 
crystalline SrFe 2 (As . 6 5Po.35)2 (Srl22P) with T c =26 K, 
which unravels the development of AFM correlations in 
the normal state and the unconventional gapless SC with 
nodes. Specific heat measurement in SC state reveals 
that the fraction of RDOS at Ep are mostly induced by 
the external field. Although the low-lying quasiparticle 
excitations are gapless in a small gap even when if =0, the 
fully-gapped s±-wave state in other bands is not affected 
significantly. Such multiband effect plays a significant 
role in keeping T c relatively high in this compound. 

Single crystalline samples of SrFe2(Asi^ 2; P a ;)2 (£=0.35 
and x=0.5) were synthesized by self-flux method[7|, |20j. 
We used small pieces of single crystals that are supercon- 
ducting at T c =26 K for £=0.35 with a sharp transition 
width AT C <1 K. 31 P-NMR measurement for £=0.35 was 
performed at fields f/=l— 14.0 T perpendicular to the 
c-axis using the aligned single crystalline samples. The 
Knight shift 31 K was measured with respect to a reso- 
nance field in H3PO4 ( 31 K=Q in Fig.[T]). For comparison, 
these measurements were also performed on £=0.5 with 
T c =8 K in the overdoped regime. 

Figure [1] shows a temperature (T) dependence of 31 P- 
NMR spectrum at iJ=14.0088 T for £=0.35. The NMR 
spectral shape does not change down to T C (7J)=23 K 
with a full-width at half maximum (FWHM) of 55 kHz 
at -14 T and 250 K, which corresponds to AK < 0.02%, 
indicating that the present single crystals are totally ho- 
mogeneous. Upon cooling below T c , its peak shifts to a 
low frequency side and the FWHM increases owing to the 
appearance of SC diamagnetism||8]. As indicated in the 
inset of the figure, 31 K is nearly T-invariant in the nor- 
mal state, but decreases below T C (H). In the £=0.5, the 
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FIG. 1: (Color online) (a) T dependence of 31 P-NMR spec- 
tra for Srl22P(:r=0.35) at #=14.0088 T. (b) T dependence 
of Knight shift 31 K, which is T-invariant, but decreases be- 
31 K for £=0.5 is slightly smaller than 



low T C (H)=23 K. The 
that of £=0.35 in the normal state, 
subtle reduction of DOS at J5f 12111 



which is attributed to the 



FWHM becomes 95 kHz at 100 K, being broader than 
that of £=0.35, due to the higher doping level. The 31 K 
in the normal state is smaller than that of a;=0.35, which 
is attributed to the subtle reduction of DOS at Ep , being 
in good agreement with the result in Bal22P[2lj|. 
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FIG. 2: (Color online) Nuclear spin-lattice relaxation rate 
1/Ti was obtained from the recovery of nuclear magnetization 
by fitting with a simple exponential recovery curve of m(t) = 
[M -M(t)]/M = exp(-t/Ti) in (a) normal and (b) SC state 
for £=0.35 at H ~14 T. 

Nuclear spin-lattice relaxation rate 1 jT\ for both sam- 
ples was obtained from the recovery of nuclear magnetiza- 
tion by fitting with a simple exponential recovery curve of 
m(t) = [M - M(t)]/M = cxp(-i/Ti) for 31 P{I = 1/2). 
Here Mo and M(t) are the respective nuclear magneti- 
zations of 31 P for the thermal equilibrium condition and 
at a time t after the saturation pulse. Typical recovery 
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FIG. 3: (Color online) T dependence of 31 P-NMR (TiT)" 1 in 
the normal state for Srl22P along with the results for Bal22P 
with T c =30 kJHJ and for BaK122 with T c =38 K |J. Here, 
the data for 31 P( 57 Fe)-NMR (TiT) -1 in Bal22P(BaK122) are 
normalized by the value at 250 K. Solid curve is a fit for 
Srl22P(£=0.35) on the assumption of (TiT) -1 = a/(T+6) + b 
with parameters a=8.7, 6=0.38, and r9= — 0.5, indicating the 
closeness to the AFM QCP. 



curves of x=0.35 are presented for normal- and SC state 
in Figs. Ufa) and^b), respectively. The m(t) can be 
determined by a single component of T\ for both s=0.35 
and 0.5. Figure [3] shows the T dependence of 31 P-NMR 
1/T\T in the normal state, which increases significantly 
upon cooling down to T c . In general, 1/T\T is described 
as, 
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where A q is a wave-vector (q)-dependent hypcrfine- 
coupling constant, xil^) a dynamical spin suscepti- 
bility, and wo an NMR frequency. Since the Knight 
shift does not change in the normal state, the increase 
of 1/TiT upon cooling is due to the development of 
AFM correlations as well as the case in Bal22P[21] and 
BaK1220. We assume two-dimensional (2D) AFM cor- 
relations model that predicts a relation of \fT\T oc 
Xq (T) oc 1/(T + 6) when a system is close to an AFM 
quantum critical point (QCP) 22]. Here, the staggered 
susceptibility xq(T) with an AFM propagation vector 
q—Q follows a Curie- Weiss law. Since 1/TiT diverges 
towards T — s> when 6 = 0, 9 is a measure of how close 
a system is to an AFM QCP. Actually, as shown by the 
solid curve in Fig. [3l the T dependence of \jT\T for 
x=0.35 can be fitted by assuming l/T x T = a/ {T + 6) +b 
with parameters a— 8.7, 6=0.38, and 9=— 0.5K. In par- 
ticular, it is notable that is nearly zero for £=0.35, 
indicating the closeness to the AFM QCP as well as in 
Bal22P with a -7, 9 - K, and b -0.17 Such a 
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TABLE I: Superconducting characteristics and lattice param- 
eters for SrFe 2 (As .65Po.35)2 (Srl22P)[3 along with those 
of Ba .6Ko.4Fe 2 As2 (BaK122) Q and BaFe 2 (As . 6 7Po.33)2 
(Bal22P)0,||. 
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Type of SC Gap 


Full 


gapless 


gapless 


o[A] 


3.917 


3.92 


3.92t 


c[A] 


13.3 


12.8 


12.23t 


hp n [K] 


1.38 


1.32 


1.32t 



f) The values for £=0.35 are interpolated from the data in refill 



critical enhancement of \/T\T disappears in x=0.5, as 
shown in Fig. [3J 

We address why the enhancement of 1/TiT at low 
T becomes distinct in going from Srl22P, Bal22P, to 
BaK122 as seen in Fig. [3] As listed in table HI the iso- 
valent substitution of P for As largely reduces the c-axis 
length and a pnictgen height (hp n ), whereas the a-axis 
length remains constant, which may be equivalent to ap- 
plying an uniaxial stress along the c-axis for (FePn) - 
layers. The band-structure calculation for Bal22P has 
predicted that one of the hole FSs mainly derived from 
the XZ/YZ/Z 2 orbitals exhibits three dimensional (3D) 
features when hp n becomes lower in the 122 compounds, 
whereas the hole FS from X 2 — Y 2 orbital keeps a quasi- 
2D character (ljj. Although the local lattice parameters 
of (FePn)- resemble for both Bal22P and Srl22P, the 
reduction of c-axis length is more significant for Srl22P. 
In this context, appearance of three-dimensionality in 
the electronic structure may cause the nesting of FSs to 
weaken in some of the multiple bands. 

Next, we present the SC characteristics in 
Srl22P(a;=0.35) with T c =26 K. Figure Ufa) shows 
a plot of 31 P-NMR Ti(T c )/Ti versus T/T c for Srl22P 
along with the results for Bal22P[8] with T c =30 K 
and BaK1220 with T c =38 K. The 1/T\s for Srl22P 
(x=0.35) decrease steeply below T c without a coherence 
peak and follow a T-linear dependence below T/T c <0.5 
when H = 14 T. The 1/TiT normalized by the value 
at T c remains a finite, as shown in Fig. 0|b). Although 
it shows a weak H dependence, it still remains a finite 
even in low H limit at T ~ 0.2T C , as seen in the inset of 
Fig. IDJb), pointing to a gapless SC dominated by a large 
contribution of low- lying quasiparticle excitations at Ep. 
It is in contrast with the power-law behavior without 
T-linear dependence in 1/Ti in BaK122 despite of the 
application of high external field (12 T), which was 
consistently accounted for by the fully-gapped s±-wave 
stated . Generally, the 1/TiT in the gapless SC state is 
related to the square of RDOS at Ep (N 2 CS ), and hence 
the fraction of RDOS (N res /No) to a normal-state DOS 
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FIG. 4: (Color online) (a) Plots of 31 P-NMR Ti(T c )/Ti and 
(b) (TiT)- 1 /(TiT)^ 1 versus T/T c for Srl22P along with those 
for Bal22P with T C C =30 K @] and BaK122 with T c =38 K 
[3- Inset shows the H dependence of (TiT)" 1 /(TiT)y c 1 at 
T ~ 0.2T C . Dotted and Solid curves are simulations on the as- 
sumption of (c) DOS for simple nodal gap with RDOS (model 
A) and (d) DOS for multigap composed of a small nodal gap 
with RDOS and full gaps (model B). In Srl22P, the model B 
is appropriate (see the text). 



Nq is given by 
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1 CTrlo 
(TiT)yi T 



Using this relation, it is deduced that N res /No is approx- 
imately 0.72 when H = 14 T. Although the 1/TiT does 
not exhibit a constant at low temperatures when H < 2.9 
T, if we assume it by a value at T ~ 0.2T C , the N ios /NqS 
are as large as 0.56 and 0.66 in H=l T and 2.9 T, re- 
spectively, being much larger than N Tes /No ~0.34(H ^4 
T) for the nodal SC in Bal22P 0, and N rcs /N =0 for 
fully gapped SC in BaK1220. The H dependence of 
RDOS in this compound is discussed by specific heat 
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measurements later. The solid curves in Figs. SJa) and 
|4]Jb) are simulations when we assume two possible SC- 
gap models: (A) single nodal gap of A((j)) = A o sin(20) 
with N les /N (denoted as model A), which was tenta- 
tively assumed for Bal22P@|, and (B) multigap with 
the fully gapped s±-wave state in some bands [J] and 
gapless in other bands under H (model B), which re- 
produce the T dependence of 1/Xi in the gapless SC 
state. The simulations based on the respective models 
A and B for Srl22P(Bal22P@) are shown by the dot- 
ted curves with parameters of 2Ao ~ 5ksT c (6ksT c ) 
and iV res /iVo=0.72(0.34) and by the solid curves with 
parameters of 2Aq ~ 5fcsT c (6fcsT c ) for bands with 
full gaps and a smearing factor of DOS 77 = 0.05Aq 
assumed in BaK1220 besides the gapless band with 
jV res /JV =0.72(0.34). The assumed DOS for the model A 
and B are indicated in Figs.|4|c) and[4|d), respectively. 
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FIG. 5: Field dependence of specific heat C(T)/T at 1.8 K 
and the 7(= liniT^o C(T)/T) estimated from the T depen- 
dence of C(T)/T above 1.8 K[i|. The solid curves are a 
fitting with the function of H a (a ~ 0.7). The broken line is 
the C(T)/T at T for BaFe 2 (Aso.7Po.3)2[I3. 



To reveal the origin of the RDOS, we present the field 
dependence of the specific heat C(T) /T at 1.8 K and the 
residual 7(= lim T ^ C(T)/T) extrapolated to T = K 
by fitting with C(T)/T = 7 + f3T 2 + 8T A in Fig. [Spot - 
The result reveals that 7 is proportional to ~ H 7 for 
< H < 14 T, indicating that the low-energy quasi- 
particle excitation is dominated by the nodal gap. The 
significant difference from the H 5 dependence expected 
in the d-wave line node gap state suggests the multi- 
band effect with the unequal sizes of gaps (24] . in addi- 
tion to the Doppler shift of the quasiparticle energy in 
the vortex state, as seen in Bal22P [|, 0, EH, H3 ■ Using 
7„ = C(T C )/T C ~ 28 mJ/ (mol K 2 )'at T c , the fraction 
of RDOS was estimated to be ~ 60% of Nq in the ap- 
plication of H=14 T, which is roughly consistent with 



the result obtained by the 31 P-NMR experiment. Even 
when H = without vortices, the residual 7 that derives 
from the impurity scattering remains ^6 mJ/(mol K 2 ) 
at T=0 K [20(, corresponding to 7/7^ ~ 0.21, which was 
larger than 0.11 in Bal22P[lCj, |ll(. Taking into account 
that the highly homogeneous electronic state is guaran- 
teed by the narrow 31 P-NMR linewidths, the low-energy 
quasiparticle DOS in the nodal gap is sensitively induced 
even by very small amount of impurity scattering. This 
behavior is in contrast with the case of BaK122 (T c = 38 
K) characterized by a fully gaped SC state, where the ap- 
plication of field of H = 6 ~ 12 T does not induce such 
low-energy quasiparticle excitation 0, [HI . It is therefore 
likely that the presence of RDOS is inherent in the 122 
iron-pnictides in which a lattice contraction takes place 
by either the chemical substitution of P for As or the ap- 
plication of pressure. The results are consistently argued 
with the multigap feature that possess a large gap-size 
ratio A^/As of a large gap (Al) to a small gap (As). 
Here, we remind the well-established cases of the multi- 
band SC with a large gap-size ratio Af,/Ac ~ 3, which 
was demonstrated in MgB-^f^ and NbSe2[27]]. Even in 
these conventional s-wave superconductors, low- lying ex- 
citations are induced due to the derealization of quasi- 
particles within the vortices at significantly low fields at 
~ -ff C 2/10. The rapid increase in N ibs /Nq at low H in 
these compounds is a consequence of the small magni- 
tude of the second gap and the associated small conden- 
sation energy. Likewise, in Srl22P, it might be expected 
that the large fractional RDOS is promptly induced by 
fields sufficiently smaller than H c2 ~ 70 T0. It can 
be attributed to that the presence of a small gap with 
nodes causes the quasiparticles within the vortices to be 
more significantly delocalized than in conventional fully- 
gapped superconductors. In this context, it is a multi- 
gap effect with a gap-size ratio of different bands being 
significantly large to induce the large contribution for 
the gapless low-lying quasiparticle excitations in the SC 
state for Srl22P. Under these situations, we claim that 
the model B (Fig. BJd)), which deals with the multigap 
system, is appropriate for Srl22P, although both models 
seem to be able to reproduce the experimental results. 
Such multiband scenario gives a reasonable explanation 
for the fact that T c = 26 K of Srl22P is still relatively 
high irrespective of the gapless SC state in some of the 
multiple bands, which may be characteristic for P-doped 
122 compounds. 

Finally, in order to shed light on an origin of the 
small gap with nodes, we compare the SC characteris- 
tics and the lattice parameters for Srl22P, Bal22P and 
BaK122 in Table fl] The fully-gapped s±-wave state in 
BaK122 is characterized by the largest c-axis length and 
the highest hp n , whereas the gapless SC state in SiT22P 
is by the lowest hp n and the reduction of c-axis length, 
which is more significant than in others. Actually, the 
residual DOS was more prominent in Srl22P than in 
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Bal22P[l [13 ED, 113], as seen via both the NMR and 
the specific heat measurements, implying that the por- 
tion of the Fermi surfaces with RDOS becomes large in 
going from Bal22P to Srl22P. According to the band- 
structure calculation for Bal22P, the orbital character 
of one of the hole sheets is very sensitive to hp n when 
As is substituted by P in 122 compounds, and pointed 
out that 3D nodal structures appear in a largely-warped 
hole FS having a strong Z 2 /XZ/YZ orbital character 
in Bal22P 16], which was supported experimentally by 
ARPES|l7|. Since more prominent 3D feature in elec- 
tronic structure may be anticipated in Srl22P than in 
Bal22P, as corroborated by the normal-state properties 
discussed above, it is possible that the origin of a nodal- 
structure in some of the multiple bands is ascribed to 
the largely-warped Fermi surface. In particular, the im- 
portant point on the scenario is that the fully-gapped s±- 
wave state in the other quasi-2D bands keeps T c still high 
through the interband scattering between quasi-2D elec- 
tron and hole FSs within the k x -k„ plane in association 
with the multiband SC nature (161] . In other context, the 
monotonous decrease of the maximum T c from BaK122, 
Bal22P, Srl22P, to Cal22P is probably due to the de- 
crease of the total SC condensation energy for the fully- 
gapped s±-wave state due to the worse nesting properties 
in some of the multiple bands. Another possible scenario 
for a nodal gap structure is the closed nodal loops located 
at the fiat parts of the electron Fermi surface, which was 
proposed in Bal22P by Yamashita et of. 121] . To reveal 



the strong material dependence of SC-gap structure in 
Fe-based superconductors, the further systematic stud- 
ies on 122 compounds are desired to focus on the rela- 
tion between the local structure of FeP„ layer and the 
gap structure. 

In summary, the 31 P-NMR and specific heat measure- 
ments on single crystallines SrFe2(Aso.65Po.35)2 (Srl22P) 
have revealed an unconventional gapless SC with nodes, 
which is accounted for by the multiband feature with the 
gap-size ratio of different bands being significantly large. 
The normal-state property is dominated by the devel- 
opment of antiferromagnetic correlations pointing to the 
closeness to an antiferromagnetic quantum critical point 
for £=0.35. However, the comparison between the other 
122 compounds reveals the reduction of AFM correla- 
tions in Srl22P, implying the worse nesting property in 
some of multiple bands due to three-dimensional feature, 
which is suggested in Bal22P by the band-calculation 
and ARPES results 13, U | : The three-dimensionality ap- 
pears in the largely warped hole Fermi surface, whereas 
the other quasi-two dimensional electron and hole Fermi 
surfaces within the k x -k v plane is not significantly mod- 
ified. If we assume the nodal gap on the largely warped 
hole Fermi surface, this scenario gives a reasonable ex- 
planation for relatively high T c regardless of the gapless 
feature in this compound, where the interband scattering 
within the quasi 2D bands stabilizes the fully-gapped s±- 



wave state. The present results in association with the 
multiband nature will develop an insight into the strong 
material dependence of SC-gap structure in Fe-based su- 
perconductors. 
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